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Abstract 

The exclusive and semi-exclusive branching ratios of the r lepton hadronic 
decay modes ( hTv Ty h^7i u T , /i _ 7r 7r°i/ T , h~ > 27T°z/ T , h~ > 37r°i/ T , 2h~h + u T , 
2h h + > 2tt°u t , 3h 2h + v T and 3h 2h + > l7r°z/ T ) were measured 
with data from the DELPHI detector at LEP. 
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1 Introduction 

The r lepton, discovered in 1975 [1], is the only lepton which is sufficiently heavy to 
decay to final states containing hadrons. Predictions for the properties of such a heavy 
lepton have been made well in advance of its discovery [2]. The taus produce intermediate 
and final-state hadrons with lower backgrounds than most other low energy processes. 

This paper describes a measurement of the decay rates of the r lepton to the different 
hadronic final states as a function of both the charged-hadron and neutral-pion multi- 
plicities, with no particle identification performed on the charged hadrons. Samples of 
different r-decay final states have been selected using both "sequential cuts" methods 
and neural networks. These analyses were complementary, allowing cross-checks of the 
results and their uncertainties. 

The DELPHI detector and data sample are described in Section 2. The method used 
to determine the branching ratios is described in section 3. The techniques used to 
separate charged leptons from hadrons are outlined in Section 4.1. Section 4.2 describes 
the reconstruction of photons and neutral pions. The selection of e + e _ — > t + t~ events is 
outlined in Section 5 and the isolated r-decays are classified according to their charged- 
particle multiplicity in Section 6. The selection of r-decays as a function of the neutral 
pion multiplicity is described in Section 7 and the associated systematic uncertainties on 
the measured branching ratios are discussed in Section 8. Section 9 presents the results 
and conclusions are drawn in Section 10. 

DELPHI has previously published results on some of the decay modes measured here 
using the 1990 data sample [3]. This paper replaces those low-statistics results. Similar 
analyses performed by other LEP experiments can be found in [4]. 

2 The DELPHI Detector and data sample 

The DELPHI detector and its performance are described in detail in [5,6]. The com- 
ponents relevant to this analysis are summarised below. Unless specified otherwise, they 
covered the full solid angle of the barrel region used in this analysis (43° < 9 < 137°) and 
lay in a 1.2 Tesla solenoidal magnetic field parallel to the beam 1 . 

The charged-particle track reconstruction was based on four different detector com- 
ponents. The principal track reconstruction device was the Time Projection Chamber 
(TPC), a large drift chamber covering the radial region 35 cm < r < 111 cm. To en- 
hance the precision of the TPC measurement, track reconstruction was supplemented by 
a three-layer silicon Vertex Detector (VD) at radii between 6 and 12 cm, an Inner Detec- 
tor (ID) between 12 and 28 cm and the Outer Detector (OD) at radii between 197 and 
206 cm from the z-axis. The TPC also provided up to 192 ionisation measurements per 
charged particle track, useful for electron/hadron separation. It had boundary regions 
between read-out sectors every 60° in cf> which were about 1° wide and which were covered 
by the VD, ID, and OD. 

The main device for 7 and tc° reconstruction and electron/hadron separation, the High 
density Projection Chamber (HPC) lay between radii of 208 cm and 260 cm. It consisted 
of 40 layers of 3 mm thick lead interspersed with 8 mm thick layers of gas sampling 
volume, amounting to a minimum of about 18 radiation lengths. In the gas layers the 
ionising particles in a shower produced electrons which drifted in an electric field into wire 

1 In the DELPHI reference frame the origin was at the centre of the detector, coincident with the ideal interaction region. 
The z-axis was parallel to the e _ beam, the x-axis pointed horizontally towards the centre of the LEP ring and the y-axis 
was vertically upwards. The co-ordinates r,<j>,z formed a cylindrical coordinate system, while 8 was the polar angle with 
respect to the z-axis. 
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chambers. In these wire chambers the induced signal on cathode pads gave a measure- 
ment of the deposited charge with sampling granularity of 10 mrad x 2 mrad x 1.0 X 
in x 9 x r in the inner 4 radiation lengths and provided up to nine longitudinal sam- 
plings of the energy deposition in a shower. The spatial precision for the starting point 
of an electromagnetic shower was 1 mrad in and 2 mrad in <fi. The energy resolution 
was AE/E = 0.31/E - 44 © 0.027. 

The Hadron Calorimeter (HCAL) was the instrumented flux return of the magnet. 
It was longitudinally segmented into 20 layers of iron and limited streamer tubes. The 
tubes were grouped to give four longitudinal segments in the readout, with a granularity 
of 3.75° x 2.96° in (px9. Between the 18th and 19th HCAL layers and also outside 
the whole calorimeter, there were drift chambers for detecting the muons which were 
expected to penetrate the whole HCAL. The barrel muon chambers (MUB) covered the 
range | cos#| < 0.602 while most azimuthal zones in the range 0.602 < | cos^l were covered 
by forward muon chambers (MUF). 

The Ring-Imaging Cherenkov detector (RICH), although not used in this analysis, had 
an important effect on the performance of the calorimetry as it contained the majority 
of the material in the DELPHI barrel region. Lying between the TPC and OD in radius, 
it covered the complete polar angle region of this analysis. The amount of material 
for particles of perpendicular incidence was equivalent to 0.6 radiation lengths and 0.15 
nuclear interaction lengths. 

The data were collected in the years 1992 to 1995, at centre-of-mass energies 
between 89 and 93 GeV on or near to the Z resonance. It was required that the VD, TPC, 
HPC, MUB and HCAL subdetectors be fully operational. The integrated luminosity of 
the data sample was 135 pb^ 1 of which about 100 pb _1 was taken at ^Js m 91.3 GeV, 
near the maximum of the Z production cross-section. 

Selection requirements were studied on simulated event samples after a detailed simula- 
tion of the detector response [6] and reconstruction by the same program as the real data. 
Samples were simulated for the different detector conditions and centre-of-mass energies 
in every year of data taking and amounted to about 16 times the recorded luminosities. 
The Monte Carlo event generators used were: KORALZ 4.0 [7] for e + e~ — > t + t~ events; 
DYMU3 [8] for e+e" -> events; BABAMC [9] and BHWIDE [10] for e+e" -> e+e" 

events; JETSET 7.3 [11] for e + e~ — > qq events; BDK [12] for four-lepton final states; 
TWOGAM [13] for e + e~ — > e + e~qq events. The KORALZ generator incorporated the 
TAUOLA2.5 [14] package for modelling r-decays. 

3 Method 

In an initial step, r-decays were selected according to their charged-particle multiplic- 
ity from a high-purity Z^ t + t~ event sample. In decays containing only one charged 
particle, this particle can be either an electron, muon or hadron. In higher charged- 
particle multiplicity decays the initial charged particles are hadrons. 

After rejection of one-prong decays containing muons and electrons the following ex- 
clusive and semi-exclusive r decay modes have been isolated and their branching ratios 
measured: 

• Charged multiplicity one: 

h~v T , h'n ^, h~2-n°v T , h~ > 37r°z/ T ; 

• Charged multiplicity three: 

2h-h + u T , 2h-h + n°u T , 2h~h + > 2vrV; 
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• Charged multiplicity five: 
3h-2h + u T , 3h~2h + > 1tt°v t . 

where h is either a n or K meson. The charge conjugate decays were also included. 

The 7T° mesons were detected and reconstructed via the photons produced in the 
decay n° — > 77. This n° decay mode has a branching ratio of (98.798±0.032)%, the 
remainder decaying through the Dalitz process n — > / ye + e~. Most of these were also 
correctly identified with the conversion rejection algorithm, and the fraction lost was a 
contribution to the inefficiency 

In the definition of these channels the presence of neutral kaons (reconstructed or not) 
was not considered. For example the decay h~v T , included channels with one charged 
hadron and none, one or more neutral kaons. The presence of neutral kaons did not 
significantly affect the selection efficiency, but was accounted for in the analysis. 

Two complementary analyses were performed on each of the samples of charged mul- 
tiplicity one and three r-decays. One analysis was based on sequential cuts and the other 
on neural networks. The r-decays were classified as a function of the tt° multiplicity 
and the branching ratios were obtained taking into account statistical and systematic 
correlations. Only a sequential cuts analysis was performed for r-decays with charged 
multiplicity of five. 

The branching ratios were measured simultaneously with the following procedure. 
Candidate r-decays can be classified using an estimator such as the maximum output 
neuron from a neural network or the set of cuts of the sequential analysis. On real data 
all decays are assigned to the different classes, providing the total number of events in 
each class: N i>obs . On simulated data, a selection-probability matrix can be obtained, 
representing the probability for decay mode j to be classified as decay mode %. This matrix 
could be diagonal, but in fact most of the off-diagonal terms are non-zero. To obtain the 
Branching Ratios Bj, a maximum-likelihood fit can then be performed to constrain the 
predicted number, N itPre d, of decays in class i to N it0 i, s . N itPre( i is given by: 



where N T is the total number of produced r particles, which is left as a free parameter in 
the fit, €j is the efficiency for decay mode j of the r + r~ selection, iV ijbfcg is the estimated 
background in class % due to non-r + r~ events, and n c is the number of classes, synonymous 
with the number of decay modes if all decays are classified. In this analysis not all 
candidate r-decays were classified as a minimum level was required on the maximum 
output neuron of the neural network. Taking into account the track multiplicity, this 
led to three additional classes, corresponding to those decays which were unclassified. 
Having three classes instead of just one for all the unclassified modes, does not improve 
the precision on the measurement, but gives additional information on the comparison of 
topological and exclusive branching ratios. 

If we do not take into account these three extra classes, the problem is undetermined, 
since there are n c + 1 unknowns (the n c branching ratios and N T ) and only n c measure- 
ments. The inclusion of these three classes, corresponding to the events not assigned to 
any given class, does not help, because, despite having three additional measurements, 
the equations are nearly degenerate (the matrix is almost singular) and the resulting fit 
is highly unstable. We avoid the problem by setting an additional constraint that all 
the branching ratios add to 1. In many previous measurements an alternative procedure 
is proposed, which is not correct in the case of multiple branching ratios. Here N T is 




(1) 
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obtained from the selected r events, together with the expected efficiency (e TT ) and back- 
ground (b), with the expression JV T = 2 • ^ • (1 - 6). However, this expression needs to 
assume a priori the branching ratios to estimate the rr selection efficiency and neverthe- 
less also makes an implicit assumption on the sum of branching ratios when computing 
that efficiency. With the method described here, unexpected decays could affect the 
goodness of the fit through its x 2 an d in particular, with an excess in the extra classes 
mentioned above. 

4 Particle identification and detector calibration 

The detector response was studied using simulation together with test samples of 
real data where the identity and momentum of the particles was unambiguously known. 
Examples of such samples consisted of e + e~ — > e + e~ and e + e~ — » events, the 

radiative processes e + e~ — > e + e~7 and e + e~ — > /i + yU~7 and Compton events selected using 
kinematic constraints. Tau-decay test samples, which were selected taking advantage of 
the redundancy of the detector, were also used. An example is r — > h(mr°)v, (n>0), 
selected by tagging the 7r° decay in the HPC. This gave a pure sample of charged hadrons 
to test the response of the calorimetry, muon chambers, and ionisation loss in the TPC. 
The decays r — > jj,uu selected with the calorimeters checked the muon chambers response 
and the TPC ionisation loss. Various test samples were used to calibrate the response 
of the model of the detector in the simulation program and where necessary to correct 
observed discrepancies. 

Further details of electron, muon and charged-hadron separation in r-decays can be 
found in the analysis of the r leptonic branching ratios [15]. 

4.1 Charged particles 
4.1.1 Tracking 

The precision on the component of the momentum transverse to the beam direction, 
p t , obtained with the DELPHI tracking detectors was A(l/p t ) = 0.0008(GeV/c)~ 1 for 
particles with momentum close to 45 GeV/c. Calibration of the momentum measurement 
was performed with e + e~ — > fi + fi~ events. For lower momenta the masses of the K® and 
A were reconstructed. For intermediate momenta three body decays (e + e~ — > /x + /z~7 
and e + e~ — ► e + e~7) were used. In these cases, the true energy of the particles can 
be calculated to a good precision from energy and momentum conservation, using the 
accurate measurement of the particle direction only. The combination of all these methods 
gives an absolute momentum scale to a precision of 0.2% over the full momentum range. 

Some 3% of hadrons reinteract inelastically with the detector material before the TPC. 
These were reconstructed with an algorithm which was designed to find secondary rein- 
teraction vertices using the tracks from outgoing charged particles produced in nuclear 
interactions. This is described in detail in the DELPHI analysis of the r topological 
Branching Ratios [16], where the efficiency of the algorithm, as well as the amount of 
material in the detector in terms of nuclear interaction lengths, were studied. The ef- 
ficiency in the data was found to agree well with the simulation prediction while there 
was an overestimate by about 10% in the simulation of the number of nuclear interaction 
lengths before the TPC gas volume. The correction factors obtained have been applied 
via reweighting techniques. 
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4.1.2 TPC ionisation measurement 

The energy loss per unit path length due to ionisation, dE/dx, of a charged particle 
traveling through the TPC gave good separation between electrons and charged pions, 
particularly in the low momentum range. The dE/dx pull variable, ridB/dx' f° r a partic- 
ular particle hypothesis (j = e,7r,K,p) is defined as 

U j dE l dx m eas ~ dE / dx expU) ^ 

[ldE / dx a(dE/dx) ' 1 ' 

where dE/dx meas is the measured value, dE/dx (j) is the expected momentum depen- 
dent value for a hypothesis j and a(dE/dx) is the resolution of the measurement. It was 
required that there be at least 38 anode sense wires used in the measurement. The dE/dx 
was calibrated as a function of particle velocity, polar and azimuthal angle. The distri- 
butions in simulation were tuned to agree with test samples of real data. The relative 
precision obtained was 6.2%. Fig. 1 shows the distribution of ^ e dE / dx and of ^- dE / dx in an 
electron test sample selected using calorimetric cuts. Fig. 2 shows the same distributions 
for a hadron test sample selected from r-decays. 



4.1.3 Electromagnetic calorimetry 

The calibration of the HPC for the energy range from 0.5 GeV to 46 GeV used test 
samples of electrons in Compton events, both radiative and non-radiative Bhabha events, 
and electrons tagged by the TPC dE/dx measurement. Since no difference was found in 
the response for electrons or photons, 7 samples were also used for the calibration. This 
will be described in section 4.2.4. 

For electrons, the associated energy deposited in the HPC (in GeV), E ass , should be 
equal to the measured value of the momentum (in GeV/c), within experimental errors. 
For hadrons the energy should be lower than the measured momentum as hadrons typi- 
cally traverse the HPC leaving only a small fraction of their energy. Muons deposit only 
a small amount of energy in the HPC. 

The ratio of the energy deposition in the HPC to the reconstructed momentum, p, has 
a peak at unity for electrons and a distribution rising towards zero for hadrons. This is 
shown in Fig. 1 for samples of electrons and Fig. 2 for samples of hadrons. It was also 
observed [15] that the energy deposition for hadronic showers starting before or inside the 
HPC had to be downscaled by about 10% in the simulation to get good agreement with 
data. This is due to an underestimate of the nuclear interaction length of the material in 
some of the subdetectors. 

Electron rejection with high hadron selection efficiency was performed using the asso- 
ciated energy deposition in only the first four layers of the HPC (corresponding to 6X 
for perpendicular incidence) in which electrons deposited a significant amount of energy, 
while hadrons had a small interaction probability. This is shown in Figs. 1 and 2 for 
electron and hadron test samples from r-decays. 



4.1.4 Hadron calorimetry and muon identification 

The signature of a muon passing through the HCAL was that of a minimum-ionising 
particle, leaving a roughly constant signal corresponding to an energy deposition of ap- 
proximately 0.5 GeV in each of the four layers, and penetrating through into the muon 
chambers. Hadrons, on the other hand, typically deposited most or all of their energy late 
in the HPC, the superconducting coil, or the first layers of the HCAL, rarely penetrating 
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through to the muon chambers. The response of the HCAL to hadrons depended on the 
energy of the hadron and where in the detector it interacted. Studies of the HCAL re- 
sponse to muons showed good agreement between data and simulation. For hadrons the 
total energy deposited in the HCAL was simulated well. However the depth profile of the 
hadronic showers was not simulated well. This is attributed to cut-offs in the modelling 
of the tails of hadronic showers in the simulation program. These had a negligible effect 
on the total deposited energy but a significant effect on the depth profile of the shower. 
This effect was corrected for by artificially adding an extra layer hit in simulated hadronic 
showers according to the results obtained from a data sample of charged hadrons pro- 
duced from a tightly-selected sample of t~ — > p~v T decays. An additional HCAL layer 
with a very low energy deposition was added in (25.5 ± 0.5)% of hadronic r decays. This 
fraction and uncertainty were obtained from a fit of the simulation shower depth profile 
to the data test sample. The distribution after this correction is shown in Figs. 3c) and 
4c). 

A number of different HCAL quantities gave hadron-muon separation, such as the 
energy deposition in the outermost HCAL layer, or the total energy in the HCAL, E hcai . 
The total associated HCAL energy, shown in Figs. 3d) and 4d), was corrected, as a 
function of the number of modules and the amount of material crossed by the particle, 
in such a way that the response for muons became independent of the polar angle. 

The muon chambers typically had between two and five layers hit by a penetrating 
muon (of momentum greater than 2.5 GeV/c.) The response to muons was calibrated 
using dimuon events. The simulation gave the same muon identification efficiency as 
the data. Most hadrons and their resultant shower did not penetrate through to the 
muon chambers, especially the external muon chambers which lay completely outside the 
magnet yoke. However, because of the poor modelling of the tails of hadronic showers in 
the simulation program, the probability that a hadron of a given momentum would leave 
a signal in the muon chambers was higher in the data than in the simulation. This was 
studied using the same data sample of hadrons in tightly-tagged — > p~v T events and in 
three-prong r-decays with very low muon contamination. Corrections were applied to the 
simulation for both the inner and outer layers of muon chambers. These were obtained 
by adding extra muon chamber hits for hadrons penetrating deeply into the HCAL so as 
to obtain good agreement between data and simulation. The fraction of extra hits was 
obtained from a fit of the muon chamber hit distribution in simulation to that for the 
data test sample. Correlations with the corrections made to the number of HCAL layers 
hit were taken into account. Figs. 3 and 4, show the response of these detectors for muon 
and hadron test samples. 

4.2 Photons and neutral pions 

The reconstruction of photons and hence of n° mesons was based principally on the 
HPC. Electromagnetic showers were reconstructed using only the HPC information with- 
out any prior knowledge of charged particles reconstructed in the tracking subdetectors 
and predicted to enter the HPC. Cuts based on the shower profile in the HPC were ap- 
plied to photon candidates to reduce the rate of fake photons from the interactions of 
hadrons in the HPC. An algorithm [6, 17] was applied to individual HPC clusters to see if 
they were compatible with having been produced by a single n° decaying to two photons 
where the showers due to the two photons overlapped significantly. In addition, photons 
which had converted to e + e~ pairs in the detector material before the start of the HPC 
were reconstructed using track segments from the tracking subdetectors. 
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4.2.1 HPC shower reconstruction 

The HPC gave up to nine longitudinal energy samples on a shower. In each sample 
the energy deposition was measured with a granularity of 2 cm in r-0 and 3.5 mm in z. 
The shower pattern recognition proceeded as follows. All samplings in all nine layers were 
projected on to a cylindrical grid of granularity 3.4 mm x 1.6 mrad in z x 0. Neighbouring 
bins were then added together into a coarser grid of granularity 0.5° by 0.5° in 9 and <fi. A 
local maximum search was performed and contiguous areas were separated if a significant 
minimum was found between two local maxima. All bins connected together after this 
were grouped together into one cluster. A fit was performed to the cluster transverse 
profiles to estimate the position of the interacting particle, together with the direction 
vector of the shower within the HPC. After the shower reconstruction, charged-particle 
tracks reconstructed in the tracking system were extrapolated to the HPC and associated 
to a cluster if it was compatible with having been produced by that particle. To increase 
the efficiency for minimum-ionising particles, additional low-energy clusters could be 
reconstructed along the track extrapolation. 

The substructure of each individual HPC cluster with energy greater than 5 GeV was 
then studied to ascertain if it was compatible with arising from a (typically high energy) 
neutral pion where the two photons from the decay produced overlapping showers. 

The high granularity of the HPC allowed a measurement of the lateral dimensions of a 
cluster. For a cluster arising from two photons entering the HPC the angular separation 
of the two photons is about m 7r o/£' 7r o for symmetric pair production (the most difficult 
case). This is about 7 mrad for E n o = 20 GeV, similar to the granularity of the detector. 
To search for cluster substructure the energy deposition inside a cluster was plotted on the 
0— 9 plane with each depth layer of the cluster weighted, giving the greatest weights to the 
more central layers, which had the most spatial-separation power. This two-dimensional 
distribution of weighted charge deposition was then fitted to a dipole function, projected 
on to the main axis, and two Gaussian distributions fitted to the projected distribution. 
The invariant mass was then calculated using the estimated energy deposition in each 
Gaussian and the opening angle calculated from the fit. Some corrections estimated from 
simulation were made to account for detector binning effects and biases in the fitting 
procedure. The main background came from photons converting just before the HPC 
and which were missed by the photon conversion reconstruction algorithm. This could 
give rise to a fake 7r° signal or a triple peak substructure in the cluster which was not 
properly handled by the algorithm. Since the magnetic field deflected charged particles 
only in 0, this problem was mostly confined to clusters with the dipole axis lying within 
100 mrad of the line with constant 9 passing through the cluster barycentre. To optimise 
the 7T°— 7 separation with a single variable, a neural network was used which had as inputs 
the estimated 7r° mass, the fraction of energy in the most energetic of the two photons 
and the angle of the dipole axis in the cluster. The network had a single output neuron 
and was trained with a sample of isolated photons in simulated /i + /i~7 final states to 
give a target output of zero and on tightly tagged n° candidates in simulated t~ — > p~v T 
decays to give a target output of unity. 

Fig. 5 shows the invariant-mass distribution and neural network output for single- 
cluster candidate 7r°'s selected from a tightly-tagged p sample in two energy ranges (8 < 
E < 12 GeV and E > 12 GeV). This Figure also shows the same quantities for an 
isolated-7 test sample from fi + fi~j. 

The HPC reconstruction was studied using isolated photons in /i + /i~7 and e + e~7 final 
states. The probability to identify a single photon as 7r° is shown as a function of the 
reconstructed HPC cluster energy in Fig. 6a); on average it was (16.8±0.6)% on data and 



8 



(15.8 ±0.2)% in simulation. The efficiency of the algorithm was studied in tightly-tagged 
r-decays containing one charged hadron and a single energetic neutral HPC cluster with 
a combined mass compatible with that of a p. Simulation studies indicate that such 
a sample of HPC clusters constituted a 90.5% pure sample of tt° — > 77 decays. The 
probability to identify a 7r° is also shown in Fig. 6b) as a function of the reconstructed 
7T° energy; on average it was (69.7 ± 0.5)% in data and (69.1 ± 0.1)% in simulation. 

The probability for a photon to be reconstructed as two HPC clusters was found to be 
a factor 1.15±0.02 larger in the data, showing an excess of unreconstructed conversions in 
the material in front of the HPC. The simulation was corrected according to this factor, 
following the reweighting technique described in [16], and a corresponding systematic 
uncertainty was assigned. 

4.2.2 Converted photons 

Photons converting in the material before the HPC fell into two classes, depending on 
whether the conversion took place before or after the TPC sensitive volume. 

About 7% of photons interacted in the material before the TPC gas volume giving an 
e + e~ pair detected in some of the tracking chambers. These were reconstructed using 
the tracks observed in the TPC. A detailed study and description of the algorithm and 
its performance can be found in [16]. In simulation the efficiency to reconstruct a con- 
verted photon was found to be (68.1 ± 0.2)% in one-prong r-decays and (59.8 ± 0.4)% 
in three-prong T-decays. Good agreement between efficiencies in data and simulation 
was observed, while the simulation program underestimated by about 10% the material 
before the TPC in terms of radiation lengths. The photons obtained with this kinematic 
algorithm were in general measured more precisely than those observed in the HPC. 

A further 35%/sin# of photons converted in the outer wall of the TPC, the material of 
the RICH inner wall, liquid radiator, drift tube walls, mirrors, and outer walls, or in the 
OD. These constituted a problem for the HPC pattern recognition as there was a more 
limited possibility to reconstruct these conversions with the tracking system as only the 
OD lay outside this region. Such conversion pairs were split in the DELPHI magnetic 
field before interacting in the HPC to produce electromagnetic depositions. This created 
a two-fold problem for the neutral particle pattern recognition: a single photon could 
produce two showers in the HPC, one from each particle of the e + e~ pair. These were 
reconstructed as either one or two clusters by the HPC pattern recognition, depending 
on the spatial separation of the showers. Potentially, both cases could be misidentified 
as a 7T° — > 77 candidate. Thus the number of reconstructed photons was incorrect. In 
particular this splitting effect was important for conversions in the outer wall of the TPC 
or the inner regions of the RICH, far from the first sensitive plane of the HPC. 

An algorithm reconstructed these converted photons from the track segments in the 
OD. The OD consisted of five layers of streamer tubes with a high efficiency for observing 
a charged particle. An OD track element direction had a resolution in azimuthal angle 
of about 1 mrad and thus gave an unambiguous determination of the sign of the charge 
of a particle up to the beam momentum, if this particle originated at radii smaller than 
150 cm. If there were two such track elements of different sign of charge in the OD, 
unassociated to reconstructed charged particles in the TPC, an algorithm which assumed 
that both track elements were produced by an e + e~ pair from a common conversion point 
was run. If this common conversion point was compatible with the material structure 
in the TPC and the RICH and the OD track elements were compatible in polar angle, 
then this was regarded as a photon. If there were HPC clusters behind the OD track 
elements these clusters had to have energies which were compatible with the estimated 
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e + and e~ energies derived from the algorithm, in which case the clusters were ignored for 
further analysis. This algorithm was typically about 25% efficient. Studies of efficiency 
using radiative dimuon and dielectron events, showed the ratio of post-TPC conversion 
reconstruction efficiency in data compared with simulation was 0.95 ± 0.07, consistent 
with unity. 

4.2.3 Hadronic shower rejection 

The granularity of the HPC was used to remove many clusters of a non-electromagnetic 
origin, such as hadronic showers occurring in the HPC or before (in the RICH or OD). 
These have different profiles in the detector due to the difference between the nuclear in- 
teraction length and radiation length of lead, and the difference in the sampling efficiency 
for the different processes through which their energy is absorbed. To be accepted as a 
photon shower a cluster had to have both longitudinal and transverse profiles consistent 
with those expected for an electromagnetic deposition [6]. This requirement rejected 
most showers from hadronic interactions. In Fig. 7, the distributions of two quantities 
related to the cluster profile in the HPC, namely the number of layers and the fraction of 
energy deposited in the first four layers, are shown for candidate photons selected after 
this pattern recognition. A good agreement between data and simulation is observed. 
Because of the high momentum of the charged hadron and the proximity to the 7T°'s, 
features typical of r-decays, additional criteria were applied to reduce further the con- 
tamination from hadronic showers. Many hadronic showers were rejected by accepting 
only those clusters for which the reconstructed energy, E s h, was greater than 500 MeV. 
The quantity d 2 sh _ ch E sh had to be greater than 10 deg 2 GeV, where d sh _ ch was the opening 
angle between the cluster and the track extrapolation at the HPC inner surface. This 
variable tends to be strongly peaked at low values for fake showers originating from splits 
of hadronic showers in the HPC, because they were typically of low energy and close to 
the track entry point in the HPC, in contrast to those originating from a photon pro- 
duced in a r decay. The distribution of this quantity is shown in Fig. 8, showing good 
agreement between data and simulation. No hadronic rejection criteria were applied to 
HPC clusters which were identified as candidate n mesons with the single-shower 7r° 
algorithm, as such clusters benefited from a low background. 

In Fig. 9 the energy spectra for selected HPC clusters are shown for the maximum- and 
minimum-energy photon in a r decay hemisphere, for different numbers of reconstructed 
clusters in that hemisphere. The agreement between data and simulation is good in all 
cases for both the low-energy region and the high-energy region. 

The full photon reconstruction efficiency was studied in two steps. First, electron 
samples where the track had left a signal in the OD, with a small probability of having 
interacted before reaching the HPC, were used to estimate the shower reconstruction 
efficiency. Isolated 7 samples from radiative ee and fifi were used to check the shower 
profile cuts. The efficiency in the data was found to be (0.3 ± 0.2)% less than in the 
simulation. 

The production of fake photons from hadronic interactions was estimated from the 
data and simulation agreement in the distribution shown in Fig. 8, for small values of 
the variable, where the fake photons rate is comparable to that of the real photons. The 
simulation was found to reproduce correctly the data to a relative 3%. 
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4.2.4 Energy scale 

In addition to the previously measured electron samples, the HPC energy scale was 
studied using isolated photons in /i + /i~7 and e + e~7 final states and compton scattering 
events. In these three cases the direction is well defined and the particle energy can be 
inferred with a very good precision using kinematic constraints, independently from the 
energy measurement in the calorimeter. This allowed the HPC energy response to be 
calibrated as a function of energy. A precision of 0.5% or better was obtained on the 
energy scale throughout the entire energy spectrum. The measured energy resolution 
was a(E)/E = 0.31 x E~°- M © 0.027. 

4.2.5 Spatial resolution 

The efficiency to reconstruct electromagnetic showers close to charged hadron tracks 
and showers in the HPC is important in r-decays where the r-decay products are tightly 
collimated. To illustrate this, Fig. 10 shows the minimum angular distance between 
different types of HPC clusters: neutral clusters fulfilling the photon requirements, those 
failing them and those associated to a charged particle. The good agreement of data and 
simulation in the region of very small opening angles demonstrates that all these effects 
are simulated correctly. 

4.2.6 Neutral pions 

Fig 11 shows, as a function of ir° energy, the probability, in a simulated p sample 
from r-decays, for a 7r° to produce a given number of HPC or converted photons. The 
efficiency to observe one or more photons from one n° in the angular acceptance of the 
HPC is high, dropping below 85% only in the region below 3 GeV. 

Reconstructed neutral pions fell into four different categories. The first class (I) con- 
sisted of 7T° candidates identified with the single cluster algorithm described in Sec- 
tion 4.2.1. The second class (II) contained 7r° candidates reconstructed from pairs of 
photons identified as separate HPC clusters, while the third class (III) contained 7r° can- 
didates reconstructed from pairs of photons, of which at least one was a reconstructed 
converted photon. The 77 invariant-mass distributions for classes II and III of candidate 
7T° are shown in Fig. 12. Class I dominated for the high-energy region, the class II con- 
tributed significantly in the region below 10 GeV, while the class III had a rather flat 
energy dependence. 

The fourth class (IV) recuperated photons in single-prong r-decays where a photon 
was accidentally associated to a charged- hadron track. For r-decay hemispheres where 
the HPC cluster associated to the track satisfied the photon-candidate requirements in 
all other respects, and where there was an additional photon candidate, the HPC cluster 
was disassociated from the track, provided that the invariant mass m 77 of the 77 system 
was greater than 70 MeV/c 2 . Simulation studies indicated that such decays were pre- 
dominantly due to the n ± n u T decay mode. The m 77 distribution for this class of n° is 
also shown in Fig. 12, before the mass cut. 

Fig 13 shows the total identification efficiency as well as the probability to classify a tt° 
in each of the four categories discussed above as a function of the r energy for simulated 
p decays. 

It is important to note that many of the high energy showers, despite not being resolved 
as 7T°, are nevertheless most likely to come from a merged 7r°. This accounted for in the 
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analyses in such a way that, depending on other variables, a single shower not identified 
as 7T° by any of the above criteria could be considered as a n°. 

5 Selection of e + e~ — ► t + t~ events 

The selection of the e + e~ — > Z — > t + t~ event sample is identical to that used in [16]. 
Only a summary is given here. 

In the e + e~ — > Z — > t + t~ reaction at ^/s = , neglecting radiative effects, the 
t + and r~ are produced back-to-back. The r's each decay to one, three or five charged 
and one or more neutral particles in a tightly collimated jet. Thus a t + t~ event is 
characterised by two low-multiplicity jets which appear back-to-back in the laboratory 
frame. Because each r emits at least one undetectable neutrino or anti-neutrino, the full 
event energy is not observed in the detector. 

Background events have various signatures which enable them to be separated from the 
signal. For the e + e~ — > qq channel, the typical charged-particle multiplicity is about 20, 
and quark fragmentation produces less-collimated jets. The e + e~ — > e + e~ and e + e _ — > 
fj, + fj,~ processes give a 1 versus 1 charged-particle topology, no neutral electromagnetic 
showers, and contain the full event energy measured in the detector due to the absence of 
final-state neutrinos. Two-photon events tend to have low energy visible in the detector 
due to the loss of the e + e~ pair in the beam-pipe. Cosmic rays can be removed using 
cuts on the distance of closest approach to the interaction region. 

The data were passed through the photon conversion algorithm outlined in Section 4.2 
to give an improved estimate of the numbers of charged and neutral particles in an event. 
To ensure that the r products lay in the acceptance of the relevant subdetectors it was 
demanded that the thrust axis of the event lie within the polar-angle region defined by 
| cos#| < 0.732 and that there be at least one charged particle in the polar-angle region 
defined by |cos#| > 0.035. The event was split into two hemispheres, each associated 
to a candidate r decay, by a plane perpendicular to the thrust axis and passing through 
the centre of the interaction region. It was required that there be at least one charged 
particle in each hemisphere. 

Hadronic decays of the Z were suppressed by requiring that there be a maximum of 
eight charged particles in an event. Background from four-fermion events was reduced, 
together with a further suppression of Z hadronic decays, by requiring that the event 
isolation angle be greater than 160°. The isolation angle was defined as the minimum 
angle between any pair of charged particles which were associated to opposite r-decay 
hemispheres. Backgrounds from fj, + fi~ and e + e~ final states and cosmic rays were reduced 
by requiring that the isolation angle be less than 179.5° for events with only two charged 
particles. 

The fi + fi~ and e + e~ contamination was reduced further by requiring that both 
p ra d = \ I ~rlr + ^rf and E rad = * / + -zfc be less than unity. The variables v\ and 

v Pi p 2 V 1 2 

p 2 are the momenta of the highest-momentum charged particles in hemispheres 1 and 2 
respectively. The quantity p[ was obtained from the formula p[ = ^/ssm9 2 /(sm9i + 
sin ^ 2 + I sin(#i + 2 )\), and p' 2 by analogy with the indices 1 and 2 interchanged. The 
angles 9i and 9 2 are the polar angles of the highest-momentum charged particle in hemi- 
spheres 1 and 2 respectively. The variables E\ and E 2 are the total electromagnetic 
energies deposited in cones of half-angle 30° about the momentum vectors p\ and p 2 re- 
spectively, while Ej = cp'j, for j = 1,2. Much of the remaining background from the 
dileptonic channels came from events containing hard radiation lying far from the beam. 
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Source of 




Background 


selection 




O.lliO.Ol 


e^e 


0.40±0.07 


qq 


0.29±0.01 


e + e~e + e~ 


0.27±0.03 




0.10±0.01 


e + e~r + T~ 


0.27±0.03 


e + e~qq 


0.02±0.01 


cosmic rays 


0.05±0.01 



Table 1: Selected non-r + r backgrounds, in percent, in the total sample. 

These events should lie in a plane. Where two charged particles and a photon were visible 
in the detector, such events were removed when the sum of the angles between the three 
particles was greater than 359.8°. 

Further reduction of the four-fermion contamination was achieved by requiring that 
there be a minimum visible energy of 0.09 x in the event. Energy deposits recorded by 
the luminometers (the SAT or STIC) at angles of less than 12° from the beam axis were 
excluded from this quantity. For events with only two charged particles, the additional 
condition that the vectorial sum of the components of the charged-particle momenta 
transverse to the beam be greater than 0.4 GeV/c was applied. Two-photon events 
typically have very low values of total transverse momentum compared with r + r~ events. 

Most cosmic rays were removed by the cut on isolation angle. Further rejection was 
carried out by requiring that at least one charged particle in the event have a perigee 
with respect to the interaction region of less than 0.3 cm in the r-0 plane and that both 
event hemispheres have a charged particle whose perigee point lay within 4.5 cm of the 
interaction region in z and 1.5 cm in r-0. 

In a final step, a neural network was used to reduce the background from hadronic Z 
decays [16]. 

The efficiency of the selection was estimated from simulation to be (51.74 ± 0.04)%. 
Within the angular acceptance it was about 85%. A total of 80337 candidate e + e" — > 
t + t~ events was selected. 

The background levels were estimated from the data themselves by fitting a nor- 
malisation factor to the background contribution in variables sensitive to a particular 
background, assuming that the shape of the background was that given by simulation, 
and where possible using particle identification to isolate particular backgrounds. The 
total background was estimated to be (1.51 ± 0.10)%. The different contributions are 
shown in Table. 1. The backgrounds from h + {i~t + t~ and t + t~t + t~ final 

states were negligible. 

6 Charged-particle multiplicity selection 

The selection of r-decays according to the charged-particle multiplicity was identical 
to that carried out for the categories 1, 3 and 5 in the DELPHI measurement [16] of the 
r topological branching ratios and only a brief description is given here. In the following 
a "good" track is defined as a track with associated hits in either the TPC or OD. The 
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VD-ID tracks include not only tracks reconstructed in the VD and ID without TPC or 
OD but also particles reconstructed from the decay products of nuclear interactions in 
the detector material. 

A one-prong r-decay was defined as a r-decay hemisphere satisfying any of the follow- 
ing criteria: 

• only one good track with at least one associated VD hit, and no other tracks with 
associated VD hits; 

• only one good track, without VD or ID hits, and one VD-ID track; 

• no good tracks, and only one VD-ID track. 

three-prong r-decays were isolated by demanding r-decay hemispheres satisfying at least 
one of the following sets of criteria: 

• three, four or five good tracks, of which either two or three had associated VD hits; 

• two good tracks with associated VD hits, plus one VD-ID track; 

• one good track with associated VD hits, plus one or two VD-ID tracks pointing 
within 3° in azimuth to a TPC sector boundary. 

Candidate five-prong r-decays were selected if they satisfied at least one of the following 
topological criteria: 

• five good tracks of which at least four had two or more associated VD hits; 

• four good tracks with associated VD hits, and one other VD-ID track. 

Additional criteria were applied in the selection of five-prong r-decays due to the large 
potential background from hadronic Z decays and misreconstructed three-prong r-decays. 
The background originating from 3/i ± > i7r°z/ T final states with a Dalitz decay was 
expected to occur at a similar level to the signal. Electron-rejection criteria based on 
E ass /P and dE/dx described in Section 4.1 reduced this background by about 70%, and 
it was further suppressed by requiring that all good tracks had a reconstructed momentum 
greater than 1 GeV/c. To reject qq events it was required that the total momentum 
of the the five-prong system be greater than 20 GeV/c. Only good tracks were included 
in the calculation of this quantity. 

Table 2 contains the efficiencies of these selection requirements for the different ex- 
clusive r-decay modes and the charged particle multiplicity selections, as obtained from 
simulation and after corrections for observed discrepancies between data and simulation 
in the rate and reconstruction efficiency of material reinteractions. 

In this analysis the quality of reconstruction of the charged particle tracks, especially 
their momentum and precision of the extrapolation to the calorimeters, was important for 
identification pourposes. Thus an additional requirement was made that candidate one- 
prong r-decays should contain a "good" track. This rejected candidate r-decays recon- 
structed with only a VD-ID track or with the inelastic nuclear interaction reconstruction 
algorithm. These have been extensively studied in [16] and the necessary corrections for 
any data/simulation discrepancies were applied, and the related uncertainties estimated. 

The data sample of r-decays contained 134421 candidate one-prong decays, 23847 
candidate three-prong decays and 112 candidate five-prong decays. 

7 Selection of (semi-) exclusive r-decay modes 

Analyses using sequential cuts and neural networks identified the different decay 
modes. In both cases, the different channel selections were applied simultaneously to 
take into account statistical and systematic correlations. 
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Table 2: Estimates of the t + t~ selection- and topology-classification efficiencies, in per- 
cent, for different exclusive decay modes, as obtained from simulation. The efficiencies 
are corrected for observed discrepancies between data and simulation in the rate and re- 
construction efficiency of material reinteractions. The quoted uncertainties are from the 
simulation statistics only. When no events are classified in a given class the Poissonian 
upper bound is taken as error. Numbers smaller than 0.005% are represented in the table 
as 0.00. 
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The following decay modes were selected using sequential cuts (where h = n or K): 
h~u T , h-7i°u T , hT > 2tt°v t , 2h~h + v T , 2h~h + > ln°v T , Zh~2h + v T and Zh-2h + -K°v T . The 
neural-network analysis was only performed for the one- and three-prong decays and 
included the following additional modes: h^2n°u T , h~ >37r°z/ r , 2h~h + n°v T and 2h~h + > 
27i°u T . It also included a measurement of the electronic and muonic branching ratios. 
Although no dedicated selection is present, we also quote the branching ratio for the 
inclusive channel h~ > l7r°z/ T , obtained by adding all the modes with at least one n°. 

In this analysis there is no explicit K° rejection or identification and the selection 
efficiencies were, to first order, independent of the presence of neutral kaons. These decays 
were therefore included in the equivalent class without K°. This was done regardless of 
the K° decay (even for the decay mode r~ — > h~K°v T — > h~ir°iT is T ) or of their interaction 
in the detector. For other mesons, the decays were classified according to the number of 
charged pions, charged kaons and neutral pions except for the decay modes containing r\ 
with subsequent decay to 77 or 7r + 7r~7 and uo with subsequent decay to 7r°7. These decay 
modes are difficult to isolate from the decay modes measured in this analysis, but are 
treated as background. Their total branching ratio was [18] (0.289 ± 0.027)%, (0.266 ± 
0.027)% in one-prong decays and (0.023±0.003)% for three-prongs. The branching ratios 
have been corrected for these backgrounds. 

7.1 Sequential-Cuts Analysis 

The various hadronic decay modes were selected with the cuts described below. The 
selection efficiencies and cross-talk between channels are given in Table 3 for the one- 
and three-prong modes, together with the backgrounds from non-T + r~ sources. Table 4 
contains the analogous information for the five-prong decay modes. The analysis for 
leptonic decays is described in [15]. 

7.1.1 One-prong decays 

In the selection of t~ — > hrv T decays, the separation of a single hadron from electrons 
and muons requires the use of most of the components of the DELPHI detector. The 
detector quantities used have been discussed in Section 4.1. The main background arises 
from t~ — > h~ir v T decays where the 7r° remains undetected, due to threshold effects or 
dead regions in the calorimeter. 

It was required that the charged particle have a momentum exceeding 0.05xpb eam . The 
mean energy per layer deposited in the HCAL, Eh ca i, was used to classify the charged- 
particle tracks into candidate and non-candidate minimum-ionising particles (MIP). For 
particles consistent with a MIP, E hca i < 8 GeV, a strong muon veto was applied, excluding 
all particles which were observed in the muon chambers or the outer layer of the HCAL. 
For the non-MIP region, Ehcai > 8 GeV, with less muon contamination, a muon veto was 
applied by excluding particles only if they were observed in the outer layers of the muon 
chambers. 

For electron rejection it was required that the electromagnetic energy deposited by 
the charged particle in the first four HPC layers did not exceed 350 MeV, and that 
the dE/dx did not exceed the expected signal of a pion by more than two standard 
deviations: ^ dE / dx < 2. (This dE/dx requirement was tightened for charged particles 
near to the azimuthal boundaries between HPC modules, where the HPC criterion gave 
poor rejection.) It was also required that the charged particle was either observed in the 
HCAL or deposited at least 500 MeV in the last five layers of the HPC. 
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Table 3: For the sequential- cuts analysis, classification efficiencies, in percent, for dif- 
ferent exclusive one- and three-prong decay modes, as obtained from simulation after 
correction for the data/simulation discrepancies discussed in the text. The bottom part 
shows the backgrounds in percent in each class from non-r + r~ sources. The quoted un- 
certainties are from the simulation statistics only. When no events are classified in a 
given class the Poissonian upper bound is taken as error. Numbers smaller than 0.005% 
are represented in the table as 0.00. 
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Table 4: For the sequential cuts analysis, the top part contains estimates of classification 
efficiencies, in percent, for different exclusive five-prong decay modes, as obtained from 
simulation after correction for the data/ simulation discrepancies discussed in the text. 
The bottom part shows backgrounds from non-r + r~ sources. The quoted uncertainties 
are from the simulation statistics only. 

Hadronic r-decays containing 7r°'s were rejected by insisting that there be no candidate 
photon, reconstructed as described in Section 4.2, in a cone of half angle 18° about the 
charged particle. 

The r-decay to h~it v T was selected by requesting an isolated charged particle with an 
accompanying n° candidate. The charged particle had to have a reconstructed momentum 
greater than 0.5 GeV/c and to be incompatible with the electron hypothesis using the 
loose cut of ^E/dx < 3-5. Candidate 7r°'s were subdivided into three different classes, 
described below: 

1. two photons, where each photon was measured as a separate electromagnetic cluster 
in the HPC or was a reconstructed conversion. The photons had to be separated by 
less than 10° and the reconstructed 7r° candidate had to have a reconstructed mass 
in the range 0.04 GeV/c 2 to 0.3 GeV/c 2 ; 

2. one shower with energy greater than 6 GeV and passing the criteria described in 
Section 4.2. This may happen either when a very energetic 7r° is not recognised as 
such by the 7r° reconstruction algorithms or when one of the photons enters a dead 
region of the calorimeter or is of too low energy to be observed in the calorimeter. 
The energy of the shower was taken as the energy of the 7r°; 

3. An identified tt° as described in Section 4.2.6. 

The /i"7T° invariant-mass distribution, calculated assuming the pion mass for the charged 
particle, is shown in Fig. 14. To reduce background it was required that the reconstructed 
h~ir° invariant mass lie in the range 0.48 GeV/c 2 to 1.20 GeV/c 2 and that the angle 
between the charged-particle direction and the ir° direction be less than 20°. 

The r-decay to h~ > 2ix°v T was selected by requiring an isolated charged particle with 
two or more accompanying 7T° candidates. The charged particle had to have a recon- 
structed momentum greater than 0.5 GeV/c. The candidate 7r°'s were reconstructed as 
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described in Section 4.2.6. Furthermore, decays with only one reconstructed ir candidate 
were accepted if there was at least one well-reconstructed photon candidate (as described 
in Section 4.2) which was not used in the reconstruction of a 7r°. 

This semi-exclusive mode had little background from non-r sources or from r-decay 
modes containing electrons and muons. The background was dominated by the h~ 
decay mode. Further rejection of the background was performed by requiring that the 
invariant mass of the /i~7r°7r° system be greater than 0.8 GeV/c 2 and that the total 
reconstructed energy be greater than 10 GeV. The pion mass was assumed for the charged 
particle and the n° mass for the 7r° candidate(s). 

7.1.2 Three-prong decays 

The signature of the decay t~ — > 2h~h + v T is of three charged particles with no ac- 
companying electromagnetic showers. A candidate 2hrh + v T decay had to have three 
charged-particle tracks in a hemisphere. The vector sum of the three charged-particle 
momenta had to have a magnitude greater than 0.166 x ^fs. It was required that there 
be no reconstructed photon of energy greater than 1.5 GeV within 10° of the three- 
charged-particle system momentum vector and that the total neutral electromagnetic 
energy in a cone of half-angle 10° around the three-charged-particle system be less than 
0.3 times the momentum of the three-charged-particle system. To reject cases where a 
photon or n° was associated to a charged-particle track extrapolation in the HPC it was 
required that the total energy associated to the three tracks in the first five layers of the 
HPC be less than 0.3 times the momentum of the three-charged-particle system. 

The r-decay to 2h~h + > lir°is T was selected by requesting three charged-particle 
tracks together with a ir° candidate. The n° candidate had to lie in the barrel region, 
| cos 6*| < 0.732, within a cone of half-angle 30° about the highest-momentum charged 
particle. 

7.1.3 Five-prong decays 

The exclusive decays r~ — > 3h~2h + v T and t~ — > 3h~2h +/ K°v T were selected from the 
inclusive five-prong sample. 

Decays with a total momentum greater than 40 GeV/c, an invariant mass of the five- 
charged-particle system greater than 1.5 GeV/c 2 or in which all photons had an energy 
less than 1.5 GeV were considered as — > 3h~2h + v T decays. Otherwise the decay was 
classified as t~ — > 3h~2h + ir°is T . 

7.1.4 Results of the sequential-cuts analysis 

The branching ratios were extracted from the data with a maximum-likelihood fit as 
described in Section 3. 

The numbers of candidate r-decays in each class are given in Table 5, together with the 
branching ratio obtained. The uncertainties quoted are statistical and take into account 
correlations between different channels. 

The invariant-mass distributions of the different classes of selected decays are shown 
in Fig. 14 for all cuts applied except those directly related to the mass. 
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Table 5: For the sequential-cuts analysis, the numbers of selected events in each class and 
branching ratios obtained. The quoted uncertainties are statistical only. 

7.2 Neural-Net Analysis 

The decay modes were also selected with the help of neural networks. Two different 
neural networks were designed, one for one-prong decays and another for three-prongs. 
They are described in this section. The events were initially separated according to their 
track multiplicity and then the selection with the corresponding neural net was applied. 
For five-prongs the sequential- cut analysis described in Sec. 7. 1.3 was applied. 

7.2.1 One-prong decays 

For the one-prong decay modes, a total of 43 input variables that could help the 
identification were studied: general variables (such as neutral multiplicities, invariant 
masses, and number of identified 7r°), charged-particle variables (such as momentum, 
dE/dx, and calorimetric energies) or neutral-particle quantities (such as energy, and 
shower-profile variables). This number was reduced first using a principal-component 
analysis, removing linearly-redundant variables after testing that they did not affect the 
performance. Then, the network was trained and tested with and without variables 
which appeared to be less significant; they were removed if the results were not degraded. 
Finally, 15 variables were used as input. These variables were: 

• the total invariant mass including charged and neutral particles; 

• the number of reconstructed photons; 

• the number of reconstructed 7r°; 

• the number of reconstructed photons not linked to any 7T°; 

• the magnitude of the momentum of the charged particle; 

• the polar angle of the charged particle; 

• the azimuthal angle, modulo 15°, of the extrapolation of the charged-particle track 
to the HPC; 

• the pion hypothesis dE/dx pull variable, ^E/dx'i 

• the number of muon chamber layers with hits associated to the charged particle; 

• the number of muon chamber outer layers with hits associated to the charged par- 
ticle; 

• the total electromagnetic energy deposited in a cone of half-angle 30° around the 
charged-particle track; 

• the energy in the HPC associated to the charged particle; 

• the energy in the inner four layers of the HPC associated to the charged particle; 

• the total hadron calorimetric energy associated to the charged particle; 



20 



• the number of layers in the HCAL associated to the charged particle. 

A feed-forward neural network [19] with one input layer, one hidden layer and one 
output layer was used. The input layer had 15 neurons, each one corresponding to one of 
the variables listed above. All the input variables were normalised to the range [—1,1]. 
Several structures were tested. Finally a net with one hidden layer of 40 neurons was 
used as the optimum in terms of efficiency and simplicity. The output layer consisted of 
six neurons. The assigned target value of these neurons was +1 for the corresponding 
class and —1 for the rest. Each neuron corresponded to one of the following decay modes: 
e~z/ T z/ c ; n~v T u^, h~v T \ h~7i®v T ; h~ it® -n v T ; h~ >37r°z/ T . 

A training procedure was performed on about 3000 simulated events for each of the 
classes, optimising the network parameters to give an answer in the output layer as close 
as possible to +1 in the neuron corresponding to the generated class and —1 in all others. 

The total sample of simulated events, excluding those used for the training, was used to 
evaluate the probabilities that a given decay be identified in a given class. The selection 
efficiencies of the different classes and the misidentification probabilities are shown in 
Table 6. 

Each of the preselected one-prong decays was processed through the neural network 
and the decay was identified as belonging to the class whose corresponding output was 
largest. Events with no output neuron above zero were not classified. The number of 
events with two or more output neurons above zero was negligible. 

The distributions of the maximum value of the output neuron for each decay mode for 
all decays are shown in Figs. 15 and 16. In most cases, this shows satisfactory agreement 
between data and simulation. 

7.2.2 Three-prong decays 

The three-prong r-decay candidates selected were divided into three classes: 2h~h + v T , 
2h~h + n°v T and 2h~h + > 2n°u T . 

A simpler network was used in this case, with all the electron/muon/hadron- 
identification variables removed and the remaining variables kept, giving a total of seven 
variables: 

• the total momentum of the charged-particle system; 

• the total electromagnetic energy associated to the charged-particle tracks; 

• the total electromagnetic energy deposited in a cone of half-angle 15° around the 
momentum vector of the charged-particle system including that associated to the 
charged-particle tracks; 

• the number of reconstructed photons; 

• the number of reconstructed tt°; 

• the number of reconstructed photons not used in a reconstructed 7r°; 

• the total invariant mass. 

The photons and 7r° had to lie in a cone of half-angle 30° about the highest-momentum 
charged particle. The hidden layer had 15 neurons and three output neurons were used. 
The network was trained with 3000 events of each of the signal classes optimising the 
network as for the one-prongs, to give outputs close to +1 in the neuron corresponding 
to the generated class and —1 in the others. Here, to reduce the background from other 
decays, the network was also trained with 3000 one-prong events that fulfilled the three- 
prong selection requirements, to give answers as close to —1 in all the output neurons. 
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Table 6: For the neural-networks analysis, the top part contains estimates of classification 
efficiencies, in percent, for different exclusive one-prong decay modes, as obtained from 
simulation after correction for the data/simulation discrepancies discussed in the text. 
The bottom part shows backgrounds from non-r + r~ sources. The quoted uncertainties 
are from the simulation statistics only. When no events are classified in a given class the 
Poissonian upper bound is taken as error. Numbers smaller than 0.005% are represented 
in the table as 0.00. 
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The event classification from the output neuron values was performed in an equivalent 
way to the one-prong case. The efficiencies and background levels for the different decay 
classes are given in Table 7. 

The distributions of the largest value of the output neurons in each decay are shown 
in Fig. 16, showing in most cases good agreement in shape between data and simulation. 

7.2.3 Results of the neural-network analysis 

As explained in Section 3, a simultaneous fit for the branching ratios was performed by 
fitting the predicted number of candidate r-decays in each class to the observed number. 
In this case, the information of the neural- net output was also used in the fit, where 
the sum over classes was extended to run over classes and bins in the neural-net output. 
Only positive values of this output were taken into account for the quoted results. The 
minimum value used in the fit was varied through the full range from —1 to without 
any variation on the branching ratio obtained, beyond that expected from statistical 
fluctuations. For the five-prong case the sequential-cuts analysis was used. The num- 
bers of selected candidate r-decays in each class are given in Table 8, together with the 
branching ratio obtained. The uncertainties quoted are statistical and take into account 
correlations between different channels. Despite the fact that the fit was not minimizing 
ax 2 , it was a maximum-likelihood fit, a x 2 is calculated to estimate of the goodness of 
the fit. Accounting only for statistical errors, = 808 for 490 d.o.f. was obtained, 
with the contribution from each channel presented in Table 8. The effect of systematic 
errors on the \ 2 wm De discussed in Section 8.10. 

The invariant-mass distributions of the different classes of selected decays are shown 
in Figs. 17 and 18. 

8 Systematics 

The systematic uncertainties due to any specific source were estimated simultaneously 
for all measured decay modes in the neural-network and sequential- cuts analyses. This 
was also the case for inclusive branching ratios, where the calculated systematic errors 
accounted for the existing cancellations between the errors of the different channels in- 
volved, leading in many cases to smaller errors. 

The systematic errors were evaluated using test samples of events as discussed in 
Sections 4.1 and 4.2. Where appropriate the relevant input variables were varied by the 
observed uncertainty and the selection and fit were repeated. The variation in the results 
was taken as an estimate of the systematic effect on the branching ratios. The effects of 
the external background and the preselection efficiency were also checked. The potential 
sources of systematic uncertainties are discussed below and summarised in Table 9. 

8.1 t + t~ selection and non-r + r _ backgrounds 

The background levels from channels other than t + t~ were varied by the uncertainties 
given in Section 5 and the fit was repeated. The observed changes on the results for the 
variation in each of the background types were added in quadrature to obtain the estimate 
of the systematic error. 

The probability of identifying a hemisphere from a background event in a given class 
was checked with the electron and muon test samples described in Section 4. The qq 
background was checked with low-multiplicity qq test samples selected applying the t + t~ 



23 



true r 


Neural network decay classification 


decay mode 




3/i ± ttV t 




Unclassified 


e~v T i> e 


0.01±0.00 


0.01±0.00 


0.00±0.00 


8.11±0.06 




0.03±0.00 


O.OOiO.OO 


O.OOiO.OO 


8.92±0.06 


-k~v t 


0.08±0.01 


0.03±0.00 


O.OOiO.OO 


12.04±0.09 


K~v T 


0.07±0.03 


0.01±0.01 


0.00±0.01 


7.03±0.27 




0.12±0.05 


0.06±0.03 


0.00±0.02 


12.05±0.44 


K-K»v T 


0.05±0.09 


0.13±0.11 


0.04±0.09 


9.56±0.90 




4.03±0.27 


1.07±0.14 


0.00±0.02 


21.99±0.57 


K Kgv T 


3.90±0.59 


1.01±0.31 


0.14±0.12 


22.50±1.28 


7T-K^K°U T 


2.24±0.44 


1.30±0.33 


0.49±0.21 


18.92±1.16 


— r\ T, r V\ 

ii 2Kgv T 


6.41±1.27 


6.47±1.27 


0.77±0.45 


27.05±2.30 


— o 

7T 7TV t 


0.41±0.01 


0.73±0.02 


0.07±0.00 


14.01±0.06 


K~-k g v t 


0.44±0.08 


0.92±0.12 


0.12±0.04 


12.84±0.41 


Ti ll K^V T 


0.22±0.10 


1.15±0.21 


0.15±0.08 


12.71±0.67 


TS— 7^0 


0.79±0.29 


l.46±0.40 


O.lliO.ll 


15.34±1.20 


- rxO 
7T TT KgU T 


(J.39±(J.13 


r A r in /i 

5.45±(J.4o 


(J.41±(J.13 


2U.98±0.83 


7X - 7^0 , . 
A 7T Ag Z/ T 


1 nn±n /i 1 

1.29±U.41 


5.(JU±(J.78 


U.15±U.14 


OO 97 11 rn 
22.3 ( ±1.5U 


7T ZlX V T 


U.Z f ±U.U1 


o on_i_n n /i 
z.z(J±U.U4 


n 7r i n no 
U. /o±U.Uz 


1 v oo_i_n 1 n 


A Z7T Z/ T 


n 99-i-n i ^ 


o QC-un /IS 


D fiO-l-D 9/1 


1 Q fiO-LI O/l 


TT" TT" ^ 7 / 

7T OTT is T 


0.10±0.03 


2.05±0.12 


1.73±0.11 


17.69±0.31 


Ott Tr Jr 7/ 

^j/l ii IS j- 


78.11±0.09 


14.10±0.08 


0.24±0.01 


4.84±0.05 


J\ TT TT £/ r 


77.79±0.64 


14.21±0.54 


0.18±0.07 


5.34±0.35 


7>-- Z^ + tt-+ 7 , 

j V i v 71 £/ r 


74.53±0.80 


15.72±0.67 


0.26±0.09 


7.17±0.48 


TT" TT" 1 TT - ^ 7 / 

Z7T 71 71 t/ T 


16.51±0.16 


69.06±0.19 


3.62±0.08 


6.99±0.11 


S7T ± 27r 7/ 

O/l ^/l t^'j- 


4.31±0.24 


59.12±0.57 


24.80±0.50 


7.22±0.30 


37T ± 37rV T 


1.63±0.36 


40.66±1.41 


46.68±1.43 


5.51±0.65 


37r~27r + z/ r 


18.47±1.76 


19.64±1.80 


0.51±0.32 


2.05±0.64 


37r-27r + 7rV r 


4.68±1.93 


30.12±4.19 


3.21±1.61 


7.72±2.44 


Source 


Non-r + r backgrounds 




O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


0.24±0.02 


e + e~ 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


0.60±0.09 


qq 


0.75±0.05 


1.89±0.12 


5.11±0.67 


1.07±0.11 


4f 


0.26±0.04 


0.22±0.05 


0.82±0.33 


0.53±0.04 



Table 7: For the neural-networks analysis, classification efficiencies, in percent, for dif- 
ferent exclusive three-prong decay modes, as obtained from simulation after correction for 
the data/simulation discrepancies discussed in the text. The last column represents the 
percentage of events not classified in any of the classes by the neural network, includ- 
ing the sequential- cuts selection of five-prong modes. The bottom part shows backgrounds 
in percent in each class from non-r + r~ sources. The quoted uncertainties are from the 
simulation statistics only. When no events are classified in a given class the Poissonian 
upper bound is taken as error. Numbers smaller than 0.005% are represented in the table 
as 0.00. 



24 



decay mode 


AT T~ 

IN umber 


branching ratio 


X 2 (bins) 


e~v T v c 


25529 


17.803±0.108 


54.9 


V 55 ) 




25860 


17.350±0.104 


160.1 


(55) 


I, — 

a v T 


19212 


12.780±0.120 


68.6 


V 55 ) 


7 

n 7rV T 


34675 


26.291±0.201 


or -i 

85.1 


V 55J 


II ZTT v t 


yoU4 


o c;o/i j_n qon 


59.0 


(55) 


h->3n°u T 


1083 


1.403±0.214 


92.1 


(55) 


3h ± v T 


12176 


9.340±0.090 


152.5 


(55) 




8909 


4.545±0.106 


77.8 


(55) 


3/i ± >2vr z/ r 


1272 


0.561±0.068 


51.1 


(55) 


5h ± iy T 


96 


0.097±0.015 


0.0 


(1) 




13 


0.016±0.012 


1.7 


(1) 



unclassified 


Number 


expected 


X 2 (bins) 


1 -prong 


18558 


18857.7 


2.2 ( 1 ) 


3-prong 


1517 


1455.1 


1-6(1) 


5-prong 


3 


5.2 


1-6(1) 



Table 8: For the neural-networks analysis, numbers of selected events in each class and 
branching ratios obtained. The quoted uncertainties are statistical only. The last column 
shows the contribution of each to the total x 2 , computed with statistical errors only. In 
parenthesis is shown the number of data points used in each case. The last three lines 
compare the measured number of unclassified events with the expectation after the fit. 

selection, except the isolation cut, which was changed to 120° < 9 iso < 150°. The 
classification rates were compared between real and simulated data and the systematic 
error was estimated conservatively as the largest of the statistical error and the difference 
between both. 

8.2 Charged-particle reconstruction 

The sources of systematic uncertainty associated with the charged-particle multiplic- 
ity selection have been studied in [16]. For track reconstruction, the sources investigated 
include: the efficiencies of the different tracking subdetectors to be included on a recon- 
structed track, both for isolated tracks and for tracks in higher track density topologies; 
effects of the TPC inter-sector boundary regions; the two-track resolution of the track- 
ing system and the efficiency to reconstruct a multi-prong r-decay as a function of the 
minimum opening angle between any two particles; the candidate r charge reconstruction. 

8.3 Material reinteractions 

Uncertainties from the photon conversion reconstruction were particularly important 
for those decay modes containing 7r°'s. The effect on the branching ratios was estimated 
by varying by their uncertainties the correction factors for the reconstructed and un- 
reconstructed conversions given in Table 1 of [16], which were obtained from data test 
samples of radiative dilepton events. The resultant uncertainties are dominated by the 
contribution from the unreconstructed conversions. A similar approach was taken for the 
nuclear interactions, with the correction factors given by Table 2 in [16]. 
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8.4 Relative efficiency of exclusive modes 

Due to mass effects and decay dynamics the momentum distributions of 7r ± and K ± 
are different even for otherwise identical final states. 

To estimate the size of these effects the proportions of charged pions and kaons in a 
given decay mode were varied by an amount consistent with the uncertainties quoted in 
the Particle Data Listings [18], the selection efficiency for that class recalculated and the 
fit repeated. The change in the measured branching ratio was taken as the systematic 
uncertainty. 

Within many classes there were a number of exclusive decay modes which differ in 
K° multiplicity, and which may not have the same selection efficiency. To estimate the 
uncertainty on the measured branching ratios, the exclusive branching ratios in a given 
class were varied within the uncertainties quoted in the Particle Data Listings [18]. The 
uncertainty on the three-prong modes also included a contribution due to the decay modes 
K~-K~-n + it Q v T and K~ K + -K~it Q v r which were not included in the simulation. 

Similarly, the decays containing rj and uj mesons were varied by the uncertainties on 
the world average to obtain systematic uncertainties on the measured branching ratios. 

8.5 Decay modelling 

The uncertainties associated with the modelling of the decays involving several pions or 
kaons were estimated by correcting the efficiencies taking into account differences between 
data and simulated invariant-mass distributions. In addition, the hadronic structure of 
the 3ir final state was varied between the default TAUOLA [14] model and that obtained 
in the DELPHI analysis of the 3tt structure in r-decays [20]. For the 3ttti structure the 
parameterisation of Model 1 of [21] was used and, as a cross-check, the parameterisation 
of 37T7T used in [20] was used to reweight the distributions of the minimum opening angle. 

The charged particle(s) produced in the various r-decay modes have different momen- 
tum spectra for the different helicity states. This leads to differences in acceptance as a 
function of the r polarisation due to cuts in the t + t~ selection. This is especially the 
case for r — > ttv t , Kv t where the momentum spectra differ most between the two helicity 
states. The analysis used the result and uncertainty from the DELPHI analysis on r 
polarisation [22]. 

8.6 Trigger 

The trigger efficiency for t + t~ final states was (99.98 ± 0.01)% for events within the 
polar-angle acceptance. Studies indicated that this inefficiency was due to events where 
both r's decayed via the r — > [ivv mode. This can be extrapolated to an inefficiency of 
(6 ± 3) x 10~ 4 for the channel r~ — > fi~u T u^. The associated systematic uncertainty was 
obtained by varying the inefficiency by its error. 

8.7 Energy and momentum scale and resolution 

The HPC energy scale was altered in the simulation by the uncertainty described 
in Section 4.2.4 and the complete analysis re-performed. The changes in the obtained 
branching ratios were taken as the uncertainty. In a similar manner the simulation 
energy was smeared and the branching ratios re-estimated. This took into account, with 
the correct correlation, different effects related to the electromagnetic calorimetry: e + e~ 
rejection, t~ — > e~z/ T z/ e identification and rejection through E ass , n° identification and 
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total invariant masses. The same procedure was followed with the momentum scale and 
resolution as given in Section 4.1.1. 

8.8 HCAL, muon chambers and dE/dx 

The correction in simulation to the tails of hadronic showers in the HCAL and muon 
chambers was modified by the uncertainties derived in Section 4.1. The analysis was 
repeated, and the observed variations in the branching ratios obtained were taken as 
uncertainties. 

The fraction of extra layers added in the simulation to give better data/simulation 
agreement was varied by the uncertainty obtained in Section 4. 1 and the analysis repeated. 
The uncertainties were taken from the variations in the branching ratios obtained. For 
the tails of showers penetrating into the muon chambers, the efficiency was varied by the 
uncertainty observed in the test samples for both muons and hadrons. 

In a similar way, the response of the HCAL and muon chambers for muons was varied 
within the uncertainties obtained in 4.1 with muon test samples. 

The dE/dx was varied in simulation for each particle according to the errors in the 
tuning described in 4.1 and the analysis re-performed. The uncertainties were taken from 
the changes in the branching ratios obtained. This affected mostly the r~ — > e~u T u c and 
t~ — > h~v T branching ratios whose separation depended most on dE/dx. 

8.9 Photon and neutral-pion reconstruction 

The photon efficiency, the probability to split one photon into two, the probability to 
create fake photons from a hadron, as well as the 7r° reconstruction efficiency and fake 
probability were checked with different test samples, as described in Section 4.2. The 
different errors were propagated to the efficiency tables and the fits were repeated. The 
observed difference was taken as systematic error. 

8.10 Summary of systematic uncertainties 

A contribution to the systematic uncertainty was included for the statistical uncer- 
tainty on the components of the selection-efficiency matrices due to the finite simulation 
sample size. 

The systematic uncertainty associated with each source and for each measured decay 
mode is shown for the neural network analysis in Table 9. The total systematic error 
was calculated as the quadratic sum of these contributions, since they were essentially 
independent. The errors for the sequential- cuts analysis were similar, but slightly larger 
in general. 

An attempt to estimate the effect of systematic errors on the goodness of the fit was 
made under the following procedure. The systematic errors were estimated bin by bin 
as the observed difference in the simulated distributions when the previously discussed 
systematic effects were varied within their uncertainties. A x 2 — 397 for 490 d.o.f. was 
thus evaluated, neglecting the bin to bin correlations (which slightly underestimates the 
X 2 )- The major contributions to the x 2 reported in Table 8 came from distortions of the 
neural network output in regions far from the cut, and where the signal and background 
separation was very clear and therefore did not affect the results significantly compared 
to the quoted systematic errors. In particular, the largest contribution, arising from the 
t~ — > ji^UrU^ channel is due to the slight widening of the sharp peak on the neural 
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network output distribution (Fig. 15) caused by the HCAL and muon chamber response 
systematics, with a migration of a small fraction of events from values close to 1 to the 
region from 0.5 to I. This has a large impact on the x 2 > but very small one on the 
results, since they are nevertheless identified clearly. Similar arguments apply to the 
second largest contribution, from the r~ — > 2h~h + v T) but in this case the neural network 
peak (Fig. 15) is narrower in data. 

9 Results 

The neural-network analysis gave for all hadronic channels better precision both in 
statistics and systematics, and included more channels. Therefore the results from this 
analysis were taken as the basic measurement, while the sequential analysis (except for 
the five-prong channels) was kept only as a cross-check. However, the performance for the 
leptonic decays is slightly worse than in [15] and therefore those results are not updated. 
Taking into account the statistical and systematic correlation of the channels with one 
or several n° some inclusive branching ratios were also derived. 

The results are shown in Table 10. 

The correlation matrix for the combined statistical and systematic uncertainties is 
shown in Table 11. 

Using the world averages [18] for the channels involving K° and neglecting this con- 
tribution for channels with more than three charged pions or kaons, we can derive the 
branching ratios shown in Table 12. In this subtraction, the total error on the world 
average was added in quadrature to the systematic error of these measurements. 

The sum of the branching ratios of channels giving one-prong topologies, taking 
into account correlations and after correcting for the decay modes not included in the 
analysis,(0.266±0.027)% [18], was (85.417±0.094±0.075)%, consistent with the DELPHI 
topological one-prong branching ratio measurement [16] B x = (85.316 ± 0.093 ± 0.049)%. 
Accounting for the strong correlation (0.80) arising from the fact that the classification is 
very efficient and few events remain unclassified, these two numbers agree to 1.3 standard 
deviations. 

These results are in good agreement with the current world averages [18]. 

10 Conclusions 

The measurement of r exclusive branching ratio to final states containing up to five 
hadrons has been performed with the DELPHI detector, with identification of neutral 
pions. Different semi-exclusive branching ratios, with only a lower bound on the number 
of 7T°, were also measured for final states containing up to six hadrons. A total of 134421 
one-prong, 23847 three-prong and 112 five-prong candidate r-decays were identified. Both 
sequential-cuts methods and neural networks have been used in the selection of exclusive 
decay modes with different neutral pion multiplicity, giving compatible results. The sum 
of the one-prong exclusive modes is consistent with our previous topological measurement. 
The good agreement in the number of observed and expected events that are unclassified 
by the neural network shows no evidence of unexpected decays. 

The branching ratios obtained are summarised in Table 10. Using the world-average 
measurements for channels involving neutral kaons, this contribution was subtracted. 
The results are summarised in Table 12. 
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3- or 5-prong decay mode 
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Table 9: Contributions in units of 10 5 to the systematic uncertainties on the branching 
ratios. 
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Decay mode 


Branching Ratio(%) 


r" -> /r > 0K u u T 

T — > 11 TV ^ UA Z/ T 
t~ -> /r27r° > 0X°i/ T 

/ — > 11 ^ 17T ^_ UJa t/ r 

r- -> h~ > 2tt° > Oi^V 
7- -> /T > 3tt > 0AT z/ r 


12.780 ±0.120 ±0.103 
zo.zyi ± u.zui ± u.ioU 
9.524 ±0.320 ±0.274 

Q7 OIQ -|- n ice xn 1 1 
O 1 .Zlo IE U. 100 IE U. 1 10 

10.927 ±0.173 ±0.116 
1.403 ±0.214 ±0.224 


r" -> 3/i ± > OiCV 
r~ -> 3^77° > OiCV 
r - -> S/^Mtt > OifV 
r~ -> 3/i ± >2tt° > OifV 


9.340 ± 0.090 ± 0.079 
4.545 ±0.106 ±0.103 
5.106 ±0.083 ±0.103 
0.561 ±0.068 ±0.095 


r" -> 5/i ± > OATV 

r" -> 5/i ± >l7r° > OifV 


0.097 ±0.015 ±0.005 
0.016 ±0.012 ±0.006 



Table 10: Measured branching ratios in percent. The uncertainties are statistical followed 
by systematic. 
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Table 11: Correlation matrix of the combined statistical and systematic uncertainties. 
The last three rows show the correlation with the topological branching ratios presented 
in [16]. 
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Decay mode 


Branching Katio(%J 


r~ — > h~it v r 

— i — \ i 

r- -> /r > 27r°z/ r 
r — > ft > D7T z/ r 


(11.571 ±0.120 ±0.114) 
(25.740 ±0.201 ±0.138) 

(9.498 ±0.320 ±0.275) 
(3o.o41 ± 0.155 ± 0.127) 
(10.901 ±0.173 ±0.118) 

[lAvo ± U.Z14 ± U.ZZ4J 


r~ — > 3/l ± l' T 

r~ -> S/^MvrV 
r~ -> 3/i ± >27r°z/ r 


(9.317 ±0.090 ±0.082) 
(4.545 ±0.106 ±0.103) 
(5.106 ±0.083 ±0.103) 
(0.561 ±0.068 ±0.095) 


r~ — > hh^T 


(0.097 ±0.015 ±0.005) 
(0.016 ±0.012 ±0.006) 



Table 12: Measured branching ratios in percent after subtraction of the contributions of 
channels including K°. The uncertainties are statistical followed by systematic. 



All the results are in good agreement with the current world averages [18] and have 
similar errors to the most precise single measurements. 
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Figure 1: Distributions for electron test samples in r-decays : a) the variable ^E/dx > b) 
the variable H^E/dx > c ) ^ e variable ^jr- ; d) the energy deposited in the innermost four 
layers of the HPC. Data are shown as dots and simulation by a solid line. The shaded 
region represents the simulated background including other r-decays. 
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Figure 2: Distributions for hadron test samples in T-decays of electron-hadron separation 
variables: a) the variable R dE / dx / b) the variable R dE / dx ', c) the variable ; d) 
the energy deposited in the innermost four layers of the HPC. Data are shown as dots 
and simulation by a solid line. The shaded region represents the simulated background 
including other T-decays. 
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Figure 3: Muon-identification variables for muon test samples in r-decays: a) number of 
muon chamber hits including the outer HCAL layer; b ) number of muon chamber hits in 
the outer muon chambers; c) number of layers in the HCAL; d) corrected deposited energy 
in the HCAL. Data are shown as dots and simulation after the corrections described in 
the text, by a solid line. The shaded area represents the non-muon background. 
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Figure 4: Muon-identification variables for hadron test samples in r-decays: a) number of 
muon chamber hits including the outer HCAL layer; b ) number of muon chamber hits in 
the outer muon chambers; c) number of layers in the HCAL; d) corrected deposited energy 
in the HCAL. Data are shown as dots and simulation after the corrections described in 
the text, by a solid line. The shaded area represents the non-hadron background. 
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Figure 5: Distribution of reconstructed invariant mass in GeV/c 2 (left) and neu- 
ral-network output (right) reconstructed with the single-cluster algorithm for showers from 
the tight p sample at different energies and for showers from the isolated-^ sample from 
fi + Data are shown as dots and simulation by a solid line. 
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Figure 6: ir° -identification probability in a single shower as a function of the shower 
energy, obtained from isolated-'-/ samples (a) and it sample in tightly-tagged p decays (b). 
Data are represented by crosses and simulation by filled circles. 
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Figure 7: Distribution for candidate photon clusters for the inclusive sample of r- decays: 
Number of HPC layers hit (a); fraction of energy deposited in the four innermost layers 
(b). Data are shown as dots and simulation by a solid line. The dashed line shows in 
arbitrary normalization the distributions for showers produced from charged hadrons. 
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Figure 8: Distributions of the quantity d 2 sh _ ch E s h used for rejection of hadronic showers in 
the HPC for the inclusive r sample of photons with all identification requirements applied 
except the cut on this variable. Data are shown as dots and simulation by a solid line. 
The shaded area shows the contribution from showers induced by charged hadrons. The 
last bin shows the overflows. 
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Figure 9: Energy distributions of identified photons in the HPC for r -decays containing 1, 
2, 3 or 4 such clusters. The figures on the left represent the most energetic cluster in the 
decay and those on the right the least energetic. Data are shown as dots and simulation 
by a solid line. 
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Figure 10: Distribution of opening angle between different types of showers: a) a neutral 
shower fulfilling the 7 identification with a shower associated to a track b) two neutral 
showers fulfilling the 7 identification, c) a neutral shower failing the 7 identification with a 
shower associated to a track c) a neutral shower fulfilling the 7 identification with another 
one failing it. Data are shown as dots and simulation by a solid line. 
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Figure 11: The fraction events, as a function ofn° energy with: (top) 0,1,2, or more than 
2, photons reconstructed; (bottom) a reconstructed conversion before or after the TPC, 
for a sample of simulated non-radiative p-decays. 
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Figure 12: Distribution of 77 invariant mass for tt° candidates for two neutral showers 
(top), neutral shower-converted photon (middle) and charged-neutral shower (bottom) . 
Data are shown as dots and simulation by a solid line. 
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Figure 13: Total n° reconstruction efficiency (top) and probability to reconstruct a n° 
in any of the classes described in the text (bottom) as a function of the n° energy for 
simulated p decays. 



40 



DELPHI 




Invariant mass (GeV/c 2 ) Invariant mass (GeV/c 2 ) 





Figure 14: Invariant-mass distributions for the decays selected with sequential cuts, ex- 
cluding the cuts directly related to this variable. Data are shown as dots, simulation by a 
solid line. The shaded area shows the background prediction from simulation. 
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Figure 15: Maximum- output neuron value in one-prong analyses. For each event the 
output of the class whose output neuron is maximum is represented. Data are shown as 
dots and simulation by a solid line. The shaded area represents the background events. 



48 



DELPHI 






Figure 16: Maximum- output neuron value in one-prong and three-prong neural-net anal- 
yses. For each event the output of the class whose output neuron is maximum is repre- 
sented. Data are shown as dots and simulation by a solid line. The shaded area represents 
the background events. 



49 



DELPHI 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 

Invariant mass (GeV/c 2 ) 




Figure 17: Invariant-mass distributions for the one-prong decays selected with the neural 
network. Data are shown as dots, simulation by a solid line. The shaded area shows the 
background prediction from simulation. 
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Figure 18: Invariant-mass distributions for the three-prong decays selected with the neural 
network. Data are shown as dots, simulation by a solid line. The shaded area shows the 
background prediction from simulation. 



